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We propose a mechanism of cell motility which is based on contraction and does not require protrusion. The 
contraction driven translocation of a cell is due to internal flow of the cytoskeleton generated by molecular 
motors. Each motor contributes to the stress field and simultaneously undergoes biased random motion in the 
direction of a higher value of this stress. In this way active cross-linkers use passive actin network as a medium 
through which they interact and self-organize. The model exhibits motility initiation pattern similar to the one 
observed in experiments on keratocytes. 



Coordinated crawling-induced movements of eukaryotic 
cells involve spatial and temporal self-organization at the cy- 
toskeletal level. In particular, to achieve a motile configura- 
tion the cell must first polarize 1 1 1. While both myosin con- 
traction and actin treadmilling contribute to cell migration, 
contraction appears to be essential for polarization, moreover, 
cells may be driven by contraction only |2|. The contraction- 
dominated motility is driven by 'pullers' and can take place 
even when 'pushers' are disabled |3|. 

In this Letter we show that the positive feedback mecha- 
nism giving rise to symmetry breaking involved in contraction 
dominated motility can be interpreted as an uphill diffusion 
driven at the microscale by advection of molecular motors. 
These motors mechanically propel the actin network by in- 
flicting contraction. In turn, the network drags the motors am- 
plifying contraction and creating an autocatalytic effect |4|. 
Such coupling leads to build up of motor concentration which 
is limited by elastic stiffness, friction and diffusion, all resist- 
ing the runaway and providing a negative feedback. 

By using the term mechanotaxis we imply conceptual sim- 
ilarity of the described motility mechanism with chemotaxis. 
Each motor generates a stress field and the other motors un- 
dergo biased random motion in the direction of a higher value 
of the stress. In this way active cross-linkers use passive actin 
network as a medium through which they interact and self- 
organize. After the symmetry of the static configuration is 
spontaneously broken the resultant active motion inside the 
cell produces overall steady translocation of the cell body. 

The idea that contraction causes flow which in turn carries 
the regulators of contraction is incorporated into the hydrody- 
namic description of active fluids |5 |. In static conditions, it 
has been shown to describe peaks in concentration of stress 
activator amplified by advective influx due to active stresses 
161 . In Hi similar idea was used to describe initiation of non- 
lamellipodial motility associated with angular cortex flows. 
Heuristic models of the Keller-Segel type 1 8 1 describing po- 
larization instability in static cells with fixed length were pro- 
posed in ||9j[T0l. In most of these models, however, the ef- 
fect of contraction is obscured by the account of other mech- 
anisms, in particular, treadmilling, and the focus is on gener- 
ation of internal flow rather than on the motion of a center 
of mass. There also exists considerable literature address- 
ing spontaneous motility driven directly by protrusion |11| 
and Turing patterning 1 12 |, or studying interaction of multi- 
ple mechanisms ITSll . 

To make the physics of mechanotaxis more transparent 



we study in this Letter the simplest analytically tractable ID 
model which captures both the symmetry breaking and the in- 
duced macroscopic motion. To decouple dynamics of actin 
and myosin we assume infinite compressibility of the cy- 
toskeleton |14|. In addition to active contractility the model 
accounts for long range elastic stiffness linking the front and 
the back of the ceU |15 |. 

We show that initiation of motility in the mechanotaxis 
model is controlled by the average concentration of motor pro- 
teins. The increase of motor concentration beyond a particular 
threshold leads to a bifurcation from a static symmetric regime 
to an asymmetric traveling wave (TW) regime describing a 
moving cell. While several TW regimes may be available for 
the same value of parameters, stable TW solutions localize 
motors in the trailing edge of the cell in agreement with ob- 
servations |2|. 

The model Consider the force balance equation for a ID 
layer of an active gel in viscous contact with rigid background 
dxO' = ^v, where a{x^t) is the stress, v{x,t) is the veloc- 
ity and ^ is the friction coefficient. Following (61 [TH [161 we 
write a = rjdxV + x^, where r] is the bulk viscosity, c is the 
concentration of motors and x > is the contractile pre- 
stress (per motor). The function c(x, t) satisfies advection- 
diffusion equation dtc + dx{cv) = Ddxx^i where D is the 
diffusion coefficient. We assume that l-{t) and are the 
unknown boundaries of the cell. We also account for a mean 
field type linear elastic interaction due to membrane or cortex 
ifTSl by using the following mechanical boundary condition, 
a{l±{t),t) = -fc(L-Lo)/Lo, where = is 
the length of the cell, k is the effective elastic stiffness and 
Lq is the reference length. Since we neglect treadmilling 
we can write the kinematic boundary conditions in the form 
^± = ^(^±)- Finally, we impose zero exterior flux of motors 
dxc{l±{t)^t) = which implies that the average concentra- 
tion Co = Lq^ c(x, t)dx is conserved. 

If we now normalize length by Lq, time, by Lq/D and 
stress by k, we obtain a Keller-Segel type system 

-Zdxxcr + cr = J^c/co, 
dtC + Xdx{cdxcr) = dxxC, 

where the dimensionless constants are Z = 7^/(^Lq), % = 
k/{^D) and 'P = cox/k. If cr is expressed through the cor- 
responding Green's function, the resulting nonlocal diffusion- 
advection problem is structurally similar to the one proposed 
in 1 9 1, however the effective kernel is different. 
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The dimensionless boundary conditions for ([T]) take the 
form a{l±{t),t) = -{L{t) - 1), dxc{l±{t),t) = and 
^±(^) = ^dxcr{l±{t)^t). They imply that the motion of the 
center of the cell G{t) = {I- {t) + /+(t))/2 is governed by the 
equation: 



G{t) 



%? f^+^'^ sh{{G-x)/VZ) 



2Zco Ji_^t) sh(L/(2V^)) 



c{x^t)dx. (2) 



One can see that if the concentration distribution is symmet- 
ric then G = and the cell cannot move, which is a simple 
analog of Purcell's theorem 1 17] with spatial asymmetry re- 
placing temporal asymmetry. From ^ one can also infer that 
the maximal speed of the cell is equal to X'J'/ {2Z). In dimen- 
sional variables 1611141 this gives xLoCo/{2r]) — lOfim/min 
which is realistic f2\ . 

TW regimes. To study the traveling wave regimes we as- 
sume that both stress and myosin concentration depend on the 
moving coordinate y = x — Vt where V is the unknown cell 
velocity. We also put l± = V and L{t) = L where L is the 
unknown length of the cell. System ([T]) reduces to a single 
equation 



Zs 



s - X{L -1)=X? 



exp(s — Vy) 



Vy)dy 



(3) 



where s{y) = X [cr{y) + (L — 1)] is the unknown function. 
The presence of four boundary conditions, 5(0) = s{L) = 
and 5^(0) = s\L) = V, ensures that both parameters V and 
L can be found along with s{y). After equation ^ is solved 
the motor concentration profile can be recovered from a rela- 
tion c{y) = coexp{s{y) - Vy)/[J^ exp{s{y) - Vy)dy]. To 
simplify the description we first assume that Z = 1 1 14] which 
means that the elastic and the viscous scales in a cell are cor- 
related. We are then left with two dimensionless parameters 
X - 100 and ? - 0.1 |7, 14], where X is the measure of 
internal stiffness while IP gives the scale of motor activity. 

The initiation of motility is associated with an instability 
of a static solution of ^ with V = 0. All such solutions 
can be written in quadratures and all of them except the ho- 
mogeneous ones imply internal flow |6|. In addition to the 
regular static solutions there are also singular static solutions 
with zero length Lq = and s{y) = lim6>^o ^f{y/^)-, where 
f{u) = {X?/2)u{l - u) and u G [0,1]; moreover, for 
IP > 1/4, those are the only static configurations. Measure 
valued solutions of this type are known in related fields [TSl 
and here they describe the collapsed cells under the action of 
unbalanced contractile stresses. 

To show that motile branches with V can bifurcate only 
from homogeneous static solutions with s{y) = 0, ^ = and 



we observe that for V 

pL 



L^ = {i±VT 



40^)72, 



(4) 



equation ([3| has an integral 

L-i exp(5(^) - Vy)dy = (1 - exp(-LF))/(LF) which 

in the limit F ^ gives exp(s) = L. Since in static 
solutions s(\j) must necessarily have a constant sign, this in- 
tegral implies that s{xj) = and hence such solutions must 




Figure 1: Three families of static solutions {L^/L- and Lo pa- 
rameterized by T at IK = 2600. The bifurcation points are la- 
beled as D 1 , 1)2 , . . . , when the nontrivial bifurcated solution is motile 
{5V / 0) and by 5*1, 5*2, ... when it is static {5V = 0). Inserts show 
the eigenfunctions 6s{y/L). In the inserts solid and dashed lines dis- 
tinguish eigenfunctions with positive and negative amplitudes (6L or 
6V). 
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Figure 2: Locus of the bifurcation points in the plane. Insert 

shows a zoom on the D 1 branch around the turning point at CP = 1/4. 
The detailed bifurcation diagrams for % — 2600 and T = 0.245 are 
shown in Fig[l]and Figjsjfrom where the meaning of labels /3, 7, jS 
7 becomes clear. 



be trivial. As we show in Fig [T] there are two families of non- 
singular trivial solutions: with longer family) and shorter 
family) lengths. 
Linearization around these trivial solutions produces the 
following linear problem for the function Ss{y/L) 



Ss 



^Ss = A^By/L 



(5) 



where u = ^X^L - A = -(cj^ + P)[{2L - 
1)/{P{L - 1))SL - (1/2) W] and B = -L{uj^ + P)SV. 
Equation ^ is supplemented with four boundary conditions 
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6s{0) = 6s{l) = 0, ^5^(0) = 6s' {1) = LSV allowing one to 
find the parameters dL and SV (up to a multiplier). Problem 
([5]) has nontivial solutions if uj ^ and 



2L^(cosa; - 1) + (cj^ ^P)ujsmuj = 0. 



(6) 




Figure 3: Bifurcation diagram with IK as a parameter showing non- 
trivial solutions branching from families of homogeneous static so- 
lutions I/+ and L-. The value T = 0.245 is fixed. Solid lines show 
stable motile branches while all the dotted lines correspond to unsta- 
ble solutions. The internal configurations corresponding to branches 
indicated by numbers (1, l^ 2, 2\ etc) are shown in Fig|4| 
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Figure 4: Internal profiles associated with successive bifurcated so- 
lutions shown in Fig[3]for T = 0.245: (1,3) correspond to asymmet- 
ric motile branches while (2,4) describe symmetric static branches. 



Solutions of the characteristic equation ^ can be split in 
two families. The first family, uj = 2m7r with m a positive 
integer, corresponds to static configurations with SV = and 
6s{y/L)) = SL{1 — cos{ujy/L)). On the parameter plane 
(IP, X), see Fig|2j sub-families of bifurcational points corre- 
sponding to different m will be labeled as 51^, 52^, ... at 
constant IP where superscripts ± indicate branches L± from 
It will be also convenient to distinguish as 51* /^l**, etc. 
bifurcational points corresponding to longer/shorter static 



configurations at fixed X. The second family, defined by 
the equation tan(co'/2) = {uj /2){1 uj'^ / L^) and indexed as 
D1^^D2^^ ... corresponds to motile solutions with 5L = 0, 
Ss{y/L) = dV{sm{u{y/L-l/2))-sm{uj/2){2{y/L)-l)y, 
the notations Z)l*/Z)l** will have the same meaning as in the 
case of static solutions. The locus of the bifurcation points 
in the parameter plane (IP, X) is shown in Fig|2] Each branch 
(say, Dl) is represented by two segments (Dl^ and Dl~) 
that meet smoothly at IP = 1/4. 

To follow the bifurcated branches into the nonlinear regime 
we performed a numerical study of the equation ([3]). A bifur- 
cational diagram at fixed IP, showing both static and motile 
configurations is presented in Figj3} the corresponding in- 
ternal profiles are shown in Fig 14] We see that each of 
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Figure 5: Bifurcation diagram with IP as a parameter showing motile 
branches connecting points DV and DV* . Corresponding bifurca- 
tion points are shown in insert in Fig|2] Parameter X is fixed in each 
graph (IK = 70 and DC = 100). 



these pitchfork bifurcations gives rise to two nontrivial solu- 
tions. For instance, point D1+ is associated with two motile 
branches, the point 5*1 + - with two static branches. Each pair 
of motile solutions is symmetric with two opposite polariza- 
tion orientations corresponding to two different signs of the 
velocity. Along the first motile branch originating at 1)1+ 
motors always concentrate at the trailing edge. For the second 
motile branch originating at 1)2+ there is an additional peak in 
the concentration profile (Fig|4]). The static bifurcation point 
/S1+ gives rise to two symmetric configurations with differ- 
ent lengths and with motors concentrated either in the mid- 
dle of the cell or near the boundaries (Fig|4]). The higher or- 
der static and motile bifurcation points produce solutions with 
more complex internal patterns. For the branches bifurcating 
from the trivial configurations belonging to I/_ family, the 
picture is similar (see Fig[3]). 

In Fig|5] we show in more detail the nontrivial solutions 
originating from the motile branch Dl at two values of pa- 
rameter X corresponding to lines af3 and af3' shown in Fig [2] 
(insert). One can see that there is a single solution connecting 
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Figure 6: Internal configuration of the moving cell on the motile 
branch Dl^ showing the localization with increasing IK at T = 
0.245. 



points 1)1* and Dl** which may belong either to one family 
I/+ (aP) or to two different families I/+ and L- (a/3 ). In the 
former case the nontrivial motile branch has a turning point at 
a finite value of IP < 1/4 giving rise to a reentrant behavior 
first observed in |9 1. In this regime the increase of the average 
concentration of myosin first polarizes the cell and initiates 
motility, but then if the concentration is increased further, the 
cell get symmetrized again and stabilize in another static ho- 
mogeneous configuration. 

Non-steady transients. A study of the initial value problem 
([T]) shows that all nontrivial solutions (static and motile) are 
unstable except for the branch bifurcating at D1+. Homoge- 
neous solutions from the I/+ family and all singular static so- 
lutions from the Lq family are stable. Numerical simulations 
also suggest that as in |6, 9|, unstable multi-peaked solutions 
are long living. This behavior is reminiscent of the classical 
spinodal decomposition modeled by ID Cahn-Hilliard equa- 
tion where the coarsening process get critically slowed down 
near multiple saddle points |19l . 



In some limiting cases the mechanotaxis equations can be 
simplified but the solutions become more singular. Thus, in 
the hyperbolic limit X ^ oo (no diffusion), the number of 
nontrivial solutions grows to infinity while the solutions be- 
come measure valued. For instance, as we show in Fig|6j 
the concentration profile for the first motile branch (Dl^ 
infinitely localizes at the trailing edge. In the inviscid limit 
Z ^ the system ([T]) reduces to dtu = dxiudxu), where 
u = 1 — %yc/cQ, which is a sign-indefinite porous flow equa- 
tion exhibiting an uphill diffusion when c/cq > . If 
the cell length is fixed meaning k oo, we have IK ^ oo 
and IP ^ and it is more convenient to restore Z and use as 
a dimensionless parameter the surviving product %y which is 
proportional to the contraction-based Peclet number x/ (^0 
introduced in |(6J [TJ ; in this case all inhomogeneous solutions 
are static and can be described in quadratures. 

In conclusion, we proposed a prototypical model of a crawl- 
ing cell showing the possibility of spontaneous polarization 
leading to steady self propulsion in the conditions when con- 
traction is the only active process while treadmilling is dis- 
abled. This model complements the existing theories of po- 
larization which place emphasis on treadmilling. The model 
reduces to a Keller-Segel type system, however, here the non- 
locality is due to mechanical rather than chemical feedback. 
We obtained a variety of motile TW regimes corresponding to 
finite size self propelling active bodies with free boundaries. 
Similar to the Navier- Stokes system, where nonlocality is hid- 
den behind the incompressibility assumption, the system of 
mechanotaxis equations has quadratic nonlinearity and shows 
an infinite sequence of bifurcations as the diffusion coefficient 
goes to zero. 
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